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I study dephasing of an underdamped nanomechanical resonator subject to random mass loading 
of small particles. I present a frequency noise model which describes dephasing due to attachment 
and detachment of particles at random points and particle diffusion along the resonator. This situ- 
ation is commonly encountered in current mass measurement experiments using NEM resonators. I 
discuss the conditions which can lead to inhomogeneous broadening and fine structure in the modes 
absorption spectra. I also show that the spectra of the higher order cumulants of the (complex) vi- 
brational mode amplitude are sensitive to the parameters characterizing the frequency noise process. 
Hence, measurement of these cumulants can provide information not only about the mass but also 
about other parameters of the particles (diffusion coefficient and attachment-detachment rates.) 
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I. INTRODUCTION 

It has been demonstrated that nanoelectromechanical 
(NEM) resonators are very suitable for sensor technology. 
In particular, it has been shown that NEM resonators can 
be employed to measure mass [THS] and charge [BHH] with 
very high precision. In the case of mass measurement 
applications, NEM resonators are useful because of the 
combination of properties such as low mass density, high 
vibrational frequencies and low intrinsic losses. These 
features are present, e.g., in carbon nanotube based res- 
onators where vibrational frequencies of /o ~ 4 GHz and 
quality factors of Q ~ 10^ have been achieved [9HTT] . 

I consider an underdamped oscillator whose eigenfre- 
quency exhibits small fluctuations. The eigenfrequency 
noise, leads to (classical) dephasing of the oscilla- 

tor response. Dephasing can be observed, e.g., in a ring 
down experiment, where, for a given realization of S(i), 
the oscillations of the oscillator coordinate, q{t), are de- 
phased with respect to the oscillator response in the ab- 
sence of frequency noise. Moreover, the average oscillator 
response, {q{t))-B.{t)i decays at a rate which is larger than 
half the oscillator energy relaxation rate, 2T. 

The effect of dephasing can be also observed in the 
power spectrum or absorption spectrum of the under- 
damped oscillator as broadening and even departure from 
the Lorentzian absorption line p^l4) . This broadening 
is not a consequence of an increase of the mechanical 
losses. Studying frequency noise in nanomechanical res- 
onators is important not only to find cures for the spec- 
tral broadening in high quality factor nanoresonators |15j 
but also in the application of the latter as sensors (e.g., 
NEM-based mass measurement applications) |16|. 

The eigenfrequency noise may have different ori- 
gins. For instance, in an electrical read-out method of 
the nanobeam displacement, charge fiuctuations in the 
nanoresonator cause frequency fluctuations due to capac- 
itive coupling with gate electrodes fn^S]. Also, random 
frequency modulation of a certain vibrational mode can 
result from parametric coupling with other modes, which 
are, e.g., thermally driven [TTI ITS] . 



In the context of NEM-based mass measurement ex- 
periments, particles enter and leave the resonator at ran- 
dom points and they also diffuse along the resonator [16] . 
Since the eigenfrequency shift of a vibrational mode de- 
pends on the particle position, random mass-loading of 
particles has been usually perceived as a limiting factor 
to the precision of NEM-based mass sensors. However, 
I show that higher order cumulants of the (complex) vi- 
brational mode amplitude can be used to characterize the 
frequency noise process |14) and determine the particle 
parameters (mass, diffusion coefficient, attachment and 
detachment rates.) 

In order to quantify the effect of random mass-loading 
in the average response of nanomechanical resonators, 
I present a frequency noise model which describes the 
quasicontinuous type of frequency noise present in NEM- 
based mass sensors. This model contains, previously 
studied, frequency noise models for particle diffusion |13j 
and discrete frequency jumps |12| as limiting cases. 

Bachtold et al. have recently realized CNT-based mass 
sensors where Xe atoms and Naphthalene molecules at- 
tach, detach and diffuse along a CNT nanoresonator |16| . 
In typical NEM-based mass measurement experiments 
the incoming fiux rate of particles can be controlled but 
the initial adsorption point and particle position on the 
resonator afterwards are beyond experimental control. 

This paper is organized as follows. In section II, I de- 
scribe the frequency noise model used to calculate the 
mode susceptibility and higher order cumulants of the 
complex vibrational amplitude. In sections HI and IV, 
I discuss dephasing due to attachment, detachment and 
diffusion of a single particle and many particles, respec- 
tively. In section V, I summarize the main results of this 
paper. 



II. MODEL 

It has been demonstrated that eigenmodes of nanome- 
chanical resonators can exhibit low intrinsic losses with 
quality factors Q ~ 10^ [HlIIl]. The vibrational ampli- 
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tude, q{t), of certain eigenmode, which is weakly driven 
near resonance by a force F cos{LdFt), satisfies the follow- 
ing equation, 



p 

- 2Tq + [ujI + 2woS(<)] Q = jj cos{ujFt) 



hit) 
M 



(1) 



where ujq, M and 2r are the bare resonance frequency, 
effective mass and energy relaxation rate of the driven 
mode, respectively. Equation ([T]) also includes the eigen- 
frequency noise, S(i), and the zero- mean additive ther- 
mal noise, fxit)- It is assumed that wq S> F, A = 
- (S(i))2)^/^ \5ujI where 5uj ^ ujf - ujq is 
the frequency detuning and tc is the frequency noise cor- 
relation time. It is also assumed that the Fourier compo- 
nents of £^{t) at frequencies uj ^ 2ujq are negligibly small 
(no parametric resonance) . The relative size between the 
quantities A, F, t^^ and \Suj\ is arbitrary. 

If the above restrictions hold then we can use the 
rotating-wave-approximation (RWA) to study the dy- 
namics of q(t) at frequencies w « wq- I introduce the 
slow (complex) dynamical variables u{t) and u*{t), de- 
fined by, 

u{t) —exp{—iujpt)^iujFq{t)+q{t)^/2iujp, (2) 

and q{t) = u{t)exp{iujpt) + u* (t) exp^—iupt). After av- 
eraging out fast oscillating terms, within the oscillation 
period 2Tr/ujp, the equation of motion for u{t) is 20 



u(i) 



[r + i{Suj - E{t))]u - i 



F 



AMujp 



+ fT{t), (3) 



where fT{t) — fT{t)exp{—iLUpt)/2iLijpM. The value 
of the moments of the complex vibrational amplitude, 
{u"{t))st, do not depend on the additive thermal noise, 
/^(t), if the latter is uncorrelated with the frequency 
noise process, S(t), and the RWA holds [Eq. (3|]. 

The frequency noise becomes correlated with the addi- 
tive thermal noise process in the case that the former 
depends on the oscillator state (backaction) . For in- 
stance, backaction occurs in the diffusion-induced bista- 
bility mechanism of driven nanomechanical resonators, 
discussed in Ref. [21], and in nonlinear oscillators, where 
the vibrational frequency depends on the vibrational am- 
plitude. Below, I neglect these extra mechanisms of 
dephasing (it is assumed that the vibrational mode is 
weakly driven.) Thus, I drop frit) from Eq. ([S]). 

The frequency noise process, S(t) = S(^(t)), is as- 
sumed to be defined in terms of a Markovian stochastic 
process ^{t), which is modeled by the master equation 



(4) 



Above, p{^, t) is the probability distribution of f at time t 
and VF(^, (,') is the matrix of transition probabilities per 
unit of time between different ^-states. It is assumed that 
W does not depend on the oscillator state, u{t). Thus, 
frit) and S(t) are independent processes. 



The evolution of the joint probability distribution of 
the oscillator state and the eigenfrequency noise pro- 
cess, P, is described by the Einstein-Fokker-Planck equa- 
tion IMl, 



dtP{u,u*,^,t) = WP + du{[r + i{6uj-E{0)]uP] 



(5) 



where S(^) is a scalar function of the noise state ^. 

The stationary value of the nth moment of the complex 
vibrational amplitude u{t) is given by, 



(u"(t)),t = J2 JJ dudu* u''P,t{u,u*,0, 



(6) 



where Pstiu,u* is the stationary solution of Eq. ^ 
and x^"'H^^) is the scaled nth moment of u{t). 

The scaled nth moment, ^("^(Jw), is written as a sum 
of scaled partial moments, x^"^^; defined by, 



(7) 



x(")(M = Ex^")(e;M 



The scaled partial moments are complex quantities which 
are coupled by the equations [T3] 



l^X^"^ - "(A - i S(e))x^"^ = mx^"~^\ 



(8) 



where A = F + 1610, x^°\6uj, ^) = PstiO and PstiO is the 
stationary solution of Eq. Q . 

Similarly, the (scaled) oscillator susceptibility, x(<^'^)j 
is defined by 



X{6uj) = {4Mup/F){u*{t))st. 



(9) 



The scaled oscillator susceptibility can also be written in 
terms of scaled partial susceptibilities x(C;<5a;) [x{Soj) = 
x(Ci which are coupled by the equation that re- 
sults from complex conjugation of Eq. (|8| with n = 1. 

The limit cases of slow [tc S> F~ , A~^] and fast 
[tc <^ F~^, A~^] frequency noise have been discussed pre- 
viously [12 [T2] . Both limit cases are described by Eq. ([s]) 
and they correspond to slow (fast) diffusion and rare (fre- 
quent) adsorption and desorption events. In what fol- 
lows, I discuss the experimentally relevant case of rare 
adsorption and desorption events, and slow and fast dif- 
fusion 16J. 

The eigenfrequency noise due to random mass-loading 
in NEM resonators is quasicontinuous. It evolves in 
time either continuously, when particles diffuse along the 
nanoresonator, or discontinuously, when particles enter 
and leave the nanoresonator. Now, I describe a model 
for the frequency noise when a single particle acts on the 
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resonator. If many particles act independently on the res- 
onator then it is only necessary to model the frequency 
noise process due to a single particle (cf. appendix B). 
This frequency noise model includes the frequency noise 
models for particle diffusion |,13j and discrete frequency 
jumps [12] as limit cases. 

The state of the particle is described by two variables 
r?(i) € {0, 1} and x{t) e [-L/2,L/2], where L is the res- 
onator length. The former variable has the value of one 
(zero) when the particle is (is not) on the nanoresonator. 
The position of the particle on the nanoresonator at time 
t is x{t). The frequency noise S(t) in Eq. ([T]) is defined 
in terms of the variables x{t) and 77 (i) by the relation 



0.25 



(10) 



where (/>(a;) is the spatial profile of the driven vibra- 
tional mode, V = miOQ/2M, m is the analyte mass, 
M = J dx po(f>{x)'^ is the vibrational mode effective mass 
and po is the resonator bare mass density. 

The transition probability matrix W{x, ry; x' , rj') for the 
frequency noise process of a single particle is defined by 



dtPn=i{x,t) ^ 
dtPn=n{t) = 



-To{x)pi{x,t) + Ti[x)J{x)po(t) 
+LD{x)piix,t), 



(11) 



L/2 



L/2 



dx [ro{x)pi{x,t) - Ti{x)f{x)po{t)], 



where pi{x,t) is the probability density that the par- 
ticle is on the nanoresonator at position x and time t. 
Similarly, poit) is the probability that the particle is not 
on the nanoresonator at instant t. The incoming fiux 
distribution per particle is f{x), which is normalized as 
/ dxf{x) — 1, and Tii^o){x) is the probability per unit 
of time that the particle is adsorbed (desorbed) at the 
point X. The detachment rate follows the Arrhenius law. 
To ~ exp{—Ea/kBT), where fcsT is the thermal energy 
and Ea is the activation energy |16j . The attachment rate 
is proportional to the incoming flux of particles and to 
the cross-section of the nanoresonator. For small vibra- 
tion amplitudes, the attachment and detachment rates 
do not depend on the oscillator state, u{t). 

The operator L o [x) describes the diffusion of an over- 
damped particle, 

LDPi{x,t) = d^{U'{x)pi)/Km + Ddlpi, (12) 

where k and D are the particle friction and diffusion co- 
efficients, respectively, and U (x) is a confining potential 
acting on the particle. I also assume reflecting bound- 
ary conditions: U' {x)pi{x^t) / nm -\- DdxPi{x^t) = at 
X ~ ±L/2. In thermal equilibrium, the diffusion coeffi- 
cient is determined by the coupling between the particle 
and the nanobeam vibrational modes (phonon bath) and 
by temperature. If the excited vibrational mode is of 
low order and weakly driven, the nanobeam temperature 
and the other (not directly driven) modes are not signif- 
icantly affected by the external drive. Hence, for under- 
damped vibrational modes, it is a good approximation to 




FIG. 1: Absorption spectrum, \mx{5oj), of the fundamen- 
tal vibrational mode of a nanobeam subject to random 
mass loading from a single particle. Curves 1 (dashed) 
and 2 (solid) depict the partial susceptibilities \mxo[5u}) and 
Im J dxx\{8uj,x), and curve 3 (dotted) depicts \mx{5uj) for 
the case of a nondiffusing particle [D — 0). Inset: Curve 4 
(dotted) depicts lmx{Sio) for the limit case of an infinitely 
fast diffusing particle (telegraph noise) and curves 5 and 6 
depict ImxiSui) for diffusion coefficients D — 25roL^ and 
D = 2r oL^, respectively. Adsorption and desorption rates 
are F^ = Fq = !//30 and 2^ is the frequency shift due to a 
particle located at the mode antinode. The oscillator damp- 
ing rate is F = 0. 



consider that the diffusion coefficient independent of the 
vibrational mode state, u{t). 

In order that the RWA is valid, the noise correlation 
time, tc ~ min {LV-D, r-\r-i}, should be much larger 
than the oscillation period, 2tt/ujq. 

Note also that Eq. ([u]) can be extended to include 
other subsystems. For instance, the substrate may be 
considered as a new subsystem which exchanges parti- 
cles with the resonator subsystem through the clamping 
regions. In this case additional exchange rates Ti/g{x) 
need to be provided. 



III. DEPHASING FROM MASS LOADING OF A 
SINGLE PARTICLE 

I use the model introduced in section II for the fre- 
quency noise to calculate the oscillator susceptibility and 
also the nth moment of the complex vibrational ampli- 
tude, u{t). I consider the fundamental ffexural mode of 
a one dimensional resonator with mode shape <^(a;) = 
a/2 cos(7ra;/_L). The shift of the vibrational frequency due 
to a particle at position x is —v4>'^{x). I assume that the 
particle may be adsorbed at any point on the resonator 
with equal probability; i.e., f{x) = l/LinEq. (11). Also, 



the adsorbed particle freely diffuses along the resonator 
[[/ = in Eq. (12)] with reflecting boundary conditions, 
dxPi{x — ±L/2[t) = 0. Below, I discuss the case where 
the maximum frequency shift per particle, 2]^, is larger 
than the oscillator damping rate, F, and the adsorption 
rate, F^, and desorption rate Fq [16|. I discuss first the 
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FIG. 2: Effect of diffusion on the imaginary part of tlie third 
cumulant, ImKs, of the complex vibrational amplitude u{t). 
Curve 1 (dotted) depicts the limit case of a nondiffusing par- 
ticle {D = 0). Curves 2 (solid), 3 (dashed) and 4 (dash- 
dotted) depict the results for a particle with diffusion coef- 
ficients D = 10"^ • ToL'^,2- IQ-^ToL'^ and 4 ■ IQ-^ToL'^, re- 
spectively. Inset: Curve 5 (dotted) depicts the limit case of 
an infinitely fast diffusing particle (telegraph noise). Curves 
6 (solid) and 7 (dashed) depict the results for D = 250rol/^ 
and D = 20r oL^, respectively. Other parameters: Fi = To = 
u/30 and F = 0. 



results for the oscillator susceptibility. 

As discussed in section II, the scaled oscillator suscep- 
tibility, x((5aj), is given by 

rL/2 

xiSuj) ^ Xo{6uj) + dxxi{Suj,x), (13) 

J-L/2 

where xoi^'-^) ^^'^ Xii^'-^^^) the scaled partial sus- 
ceptibilities associated with no particle and one particle 
at position x on the resonator. The equation for the 



scaled partial susceptibilities [cf. Eq. ( 18 )] can be solved 
analytically in terms of continued fractions (cf. appendix 

Figure [l] depicts the scaled absorption spectrum, 
ImxiSoj), for the limit cases of a slow and fast diffus- 
ing particle and v S> ri,ro,r. In the slow diffusion 
limit, the absorption spectrum exhibits inhomogeneous 
broadening. For frequency detuning du! w —Si', the 
scaled partial susceptibility of a particle at position x 
is xi{Suj, x) ~ ipist{x)/[r + To — i{S(jj + i^(fP{x))] and the 
scaled oscillator susceptibility is [13] 



iSuj) (r 



i{Suj + 2i/)) 



(14) 



The scaled partial susceptibility xoi^^) contributes to 
the total oscillator susceptibility mainly for frequency de- 
tuning Slu ~ 0. This is reasonable because xo(<5w) corre- 
sponds to no adsorbed particle (zero frequency shift). As 
depicted in Fig.[l] lmxo{Suj) is approximately Lorentzian 
with a half-width F -I- F^ and height proportional the 
stationary probability of no particle on the resonator: 



FIG. 3: Sensitivity of the cumulants, ki, K2 and K3 of u{t) 
on the particle diffusion coefficient, D. The dashed black 
and solid red curves corresponds to D = S.SOuL^ and D = 
3.33uL^, respectively. The other parameters are derived from 

r+Vi = 3.333 X 10"^!/, t+Vo+u'^l^/Stt^d = 3.733 X iQ-V, 

Po,st ~ 0.505 and f is kept fixed. Note that this choice of 
parameters slightly affects the first cumulant, m. Thus, in 
order to determine all the oscillator and particle parameters, 
it is necessary to consider the higher order cumulants. 



Po,st = Fo/(Fi + To). The asymmetry observed in 
lmX(j{6ui) at Suj = is the result of coupling with the 
scaled partial susceptibilities Xi{^^ (5w « 0). 

The limit case of a fast diffusing particle is depicted 
in the inset of Fig. [ij assuming that 1/ ^ Fj,Fo,F. In 
this limit, the oscillator absorption spectrum can be ap- 
proximated by two Lorentzians centered at Slo ~ and 
6uj — —'^{4>ix)'^)pisf These Lorentzians have half-widths 
approximately equal to F -I- F^ and F -I- Fq -I- j/^L^/Stt^Z?, 
and heights approximately equal to Po,st/(r + F,j) and 
(Lpi,5t)/(F + Fo + v^L'^/Sn^D); respectively. These re- 
sults agree with dephasing due to telegraph frequency 
noise, which is obtained in the limit D 00 [13]. Thus, 
the effect of a finite but large diffusion coefficient {D ^ 
L^Fo) is to increase the half-width of the Lorentzian cen- 



tered at (5a; = —J/ by an amount 



,2t2, 



in-^D. As the 



particle diffusion gets slower, this Lorentzian gets broader 
and it eventually looses its Lorentzian shape and acquires 
the shape given by Eq. ( [l4| . 

Higher order moments of u{t) are calculated by solv- 
ing numerically the coupled system ([s]). In the charac- 
terization of the dephasing process, cumulants are more 
important than the moments of u(t) because the former 
vanish in the absence of frequency noise [H]. Figure [2] 
shows the spectrum of the imaginary part of the third 
cumulant, K3 = (u^) — 3{u) ■ {u^) + 2{u)^, for both limit 
cases of a fast and a slow diffusing particle. The mag- 
nitude and shape of the higher order cumulants spectra 
depend on the dephasing process parameters (i.e., diffu- 
sion coefficient and attachment and detachment rates); 
moreover, this dependency is more significant for higher 
order cumulants of u(t). 

Using a frequency noise model such as the one intro- 
duced in section II, measurement of the cumulants of 
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u{t) can provide enough information to characterize the 
dephasing process and determine the particle parame- 
ters. This is demonstrated in Fig. [3] where the imaginary 
parts of the second and third cumulant are shown to be 
more sensitive to a small change of the diffusion coeffi- 
cient (solid and dashed lines correspond to D ~ S.^SL^v 
and D = 3.50L^j^, respectively.) Note that the other par- 
ticle parameters and oscillator damping rate have been 
chosen such that the first cumulant is minimally affected 
(i.e., r Tj, r -h To v'^L'^/Sn'^D, v and po,.!!t are fixed 
in the solid and dashed curves of Fig. [3|. 

In order to measure the moments and cumulants 
of the slowly varying complex amplitude u(t) = 
X[t) — ^y(^)^ one has to measure the inphase, X(t\ 
and quadrature, envelopes of the oscillations: 

q{€) = 2X{t) cos{ujpt) -\- 2y(t) sm{ujpt). In carbon-based 
nanoresonators, these measurements can be performed 
using the FM mixing technique [23] . I point out that the 
recording time of X{t) and Y{t) has to be much larger 
than the correlation time of the frequency noise and the 
oscillator energy relaxation time. 



IV. DEPHASING FROM MASS LOADING OF 
MANY PARTICLES 

In this section I consider dephasing due to particles 
acting independently on the nanoresonator. Particles act 
independently on the resonator if the average particle 
density on the resonator is small such that they do not 
interact with each other. For simplicity, it is assumed 
that the particles have identical mass, adsorption and 
desorption rates and diffusion coefficient. Hence, the to- 
tal frequency noise ^(t) is equal to a sum of independent 
and identical processes £,i{t), which are realizations of the 
process (10 1, 



N 



(15) 



1=1 



where TV is the total number of particles in the system. 
Particles can be either in the gas subsystem or in the 
resonator subsystem. The number of particles on the 
resonator n{t) evolves according to a Bernoulli process 
with master equation 

Pn = -\ron + Ti{N -n)\pn + {n + l)ToPn+i + 

r,(7V-(n-l))p„_i, (16) 

where A'o = NTi/iTi -f To) is the average number of 
particles on the resonator. The stationary distribution 
of Eq. ^ for Nq <^ N and F, = {Nq/N)To < F^ is the 
Poisson distribution: p** = exp(— 7Vo)A',^Vri!. 

Figure |4] depicts the oscillator absorption spectrum for 
A'o = 3 and v ^ Fo,F. The limit of fast diffusing par- 
ticles agrees with the results obtained previously for dis- 
crete frequency jumps with a stationary Poisson distribu- 
tion. The absorption spectrum displays a fine structure 
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FIG. 4: Absorption spectrum, Imx^Su), of the fundamental 
flexural mode of a nanobeam subject to random mass loading 
of many particles. The average number of particles on the 
resonator is A'^o = 3. The dotted curve depicts the results 
for the limit case of infinitely fast diffusing particles [D — > 
oo). The solid fine depicts the result for diffusion coefficient 
D = SOOPoL^ Inset: Limit case of slowly diffusing particles 
(D 0). Other parameters: N = 10*, To = 0.002;/, F, = 
{No/N)Vo and F = 0. 



(dotted curve) and it is formed by Lorentzians centered 
at frequencies 5lo = —ju (j = 0, 1, • • • ) with half-widths 
Tj — F-f-(j + A^o)ro [{i + Nq)T o is equal to the reciprocal 
mean time of having j particles on the resonator], and 
the heights are proportional to the stationary distribu- 
tion mj. 

The effect of a finite diffusion coefficient, D, is that 
the half-width of the above Lorentzians is increased by 
an amount « jV^L^/(87r^£>), cf. Fig. [I] (sohd curve). 
In the limit of slow diffusion, I again observe inhomoge- 
neous broadening due to the randomness associated with 
the particle insertion point. This is depicted in the in- 
set of Fig. |4j Here, the absorption spectrum, Imx((5a;), 
also exhibits sharp peaks. However, these peaks are not 
Lorentzians and they get sharper as Fq decreases. They 
also appear in the case of dephasing due to a single par- 
ticle acting on the resonator, cf. Eq. (14). The above 



results agree with Monte Carlo simulations [241 [25 



V. CONCLUSIONS 

I have discussed the effect of random mass loading of 
small particles on the resonant response of the funda- 
mental flexural mode of a nanobeam. A model for the 
frequency noise has been presented which accounts for 
particle attachment and detachment at random points on 
the resonator and particle diffusion along the resonator. 
Analytical and numerical results have been presented for 
the mode susceptibility and for some higher order cu- 
mulants of the complex vibrational mode amplitude. In 
the case of rare adsorption and desorption events and fast 
diffusion, the oscillator absorption spectrum exhibits fine 
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structure. The diffusion contribution to the Unewidths of 
the (approximately) Lorentzian hues, centered at —jv, is 
« jv^L^ /{Stt'^D), where D is the diffusion coefficient, 
j = {0,1, ■ ■ ■ ), 2i/ is the frequency shift due to a particle 
located at the vibrational mode antinode and L is the 
nanobeam length. It has been also demonstrated that 
higher order cumulants are more sensitive to the param- 
eters of the frequency noise and their measurement can 
be used to develop a mass sensing scheme, which is able 
not only to determine the mass of the particle but also its 
diffusion coefficient, attachment and detachment rates to 
the nanoresonator. 



and Q{n) = T + To — i{Suj + v) + n^r^^. The correlation 
time of frequency noise due to only particle diffusion is 
TD = {L/2tt)^/D [13]. 



VII. APPENDIX B 

Here I consider the calculation of the (scaled) suscep- 
tibility function, x((5w), for the case of many particles 
acting independently on the resonator, cf. section IV. 
The scaled oscillator susceptibility can be expressed as 



VI. APPENDIX A 

Here I derive the analytical expression for the (scaled) 
oscillator susceptibility xi^'^) = {4:Mujf / F){u* (t)) for 
the case of mass loading due to a single particle, cf. sec- 
tion HI. The assumptions are uniform incoming flux dis- 
tribution f{x) ~ 1/L and free diffusion along the res- 
onator [U{x) — 0]. The equation for the scaled sus- 
ceptibility is given by the complex conjugate of Eq. rfsl) 
with n = 1 and W defined by the master equation 
whose stationary solution is po.^t = To/{Ti + To) and 
Pi,st{x) = {1/L) ■ Ti/{ri + To). The total scaled suscep- 
tibility is equal to the sum of the scaled partial suscep- 
tibilities xo{S^) (no particle on resonator) and Xi(a;; Siv) 
(particle at position x). 



L/2 



L/2 



dxxi{x;Suj), 



(17) 



and the scaled partial susceptibilities satisfy the equa- 
tions 

[A* -iiy^^{x)+To-Ddl]xi-T,f{x)xo = ipi.sti^). 

(A*+ri)xo-ro dxxi{5uj,x) = ipo,st.,{18) 



where A* = T^iSu. The solution for xi{x]Suj) is sought 
in the form 



Xi 



{x; Suj) — ^ Ak cos{2Trkx/L). 



k>0 



The coefficients satisfy a difference equation system 
which can be solved in terms of continued fractions, 



L ■ ipi.st. + i^iPo,st./{^ + T-i - i6u 



dxxi{x;Suj) 

Xo('5w) 
where R is given by 



R{D,Sio)-T,To/{T + T, 
J dxxi{6uj,x) + ipo,st. 



— iSu 



T + T, 



iSuj 



R{D,Suj) = 6(0) 



6(1) 



(19) 



(20) 



6(2) 



xiSuj) 



dte 



At „N 



it), 



(21) 



where N is the total number of particles in the system 
(in and out of the resonator), A* = F — iSco and the 
dephasing-dependent term g{t) is equal to 



g{t) = / dx^i{x,t)+^oit). 



(22) 



The time-dependent partial susceptibilities xi(a;,f) and 
Xo(^) are related to Xi{x;6uj) and xo{^^) by Fourier 
transform 



Xi{x;Suj) 



dt 



Xi{x,t), 



(23) 



and a similar relation between xo(i) and xo('Jw) holds. 
These new (scaled) partial susceptibilities satisfy the 
equations, 

dtXi = [To-tiy(f>{x)^~Ddl]ii-TJ{x)io, 
dtXo = r,^o-r„ / dxMx,t), (24) 



with initial conditions xi{xi^) — '^Pist{x) = i{l/L) ■ 
^i/{Ti + To) and Xo(0) = ipost = i^o/i'i'i + ^o) with 
reflecting boundary conditions: dxXi{x — ^L/2,t) = 0. 
Note that g{t) can be also obtained by taking inverse 
Fourier transform of the result ( 19 1 with F = 0. 
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